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Failure of amorphous solids is fundamental to various phenomena, including landslides and earth¬ 
quakes. Recent experiments indicate that highly plastic regions form elongated structures that are 
especially apparent near the maximal shear stress E max where failure occurs. This observation 
suggested that E max acts as a critical point where the length scale of those structures diverges, 
possibly causing macroscopic transient shear bands. Here we argue instead that the entire solid 
phase (E < E max ) is critical, that plasticity always involves system-spanning events, and that their 
magnitude diverges at E max independently of the presence of shear bands. We relate the statistics 
and fractal properties of these rearrangements to an exponent 6 that captures the stability of the 
material, which is observed to vary continuously with stress, and we confirm our predictions in 
elastoplastic models. 


Amorphous solids, such as emulsions, sand or molecu¬ 
lar glasses are yield stress materials: they behave as solids 
if the applied shear stress E is low, but flow as fluids if 
it is large. Unlike the melting transition, the associated 
phase transition is dynamical: the solid phase is an ar¬ 
rested, glassy state whose properties depend on prepara¬ 
tion. One such property is failure [I], which occurs as E 
increases toward a history-dependent stress X max where 
macroscopic flow starts. For densely prepared materi¬ 
als, the stress overshoots: E max > s c HE], where X c is 
the minimum stress at which flow can be maintained in 
stationary conditions. Flow then tends to localize along 
transient (but sometimes long-lasting) shear bands [4]. 
By contrast, for loosely prepared materials E max = E c [2] 
and shear-banding may be avoided |5]. Despite its impor¬ 
tance in human applications and geophysical phenomena, 
including landslides and earthquakes [6], the microscopic 
mechanisms controlling plasticity and failure remain de¬ 
bated. 

In granular materials, recent experiments [TIE] and nu¬ 
merics m support that for E < E max , plasticity occurs 
via localized rearrangements, or shear transformations 
m, which tend to organize into elongated structures 
whose magnitude grows as X —)► X max . In m it was 
argued that for a dense initial state (X max > X c ), X max 
acts as a critical point where a correlation length £ di¬ 
verges and avalanches become system-spanning, which 
may in turn trigger macroscopic shear bands. This view¬ 
point complements the growing consensus that the re¬ 
verse transition, occurring when flows stop as the stress 
is decreased toward E c , is accompanied by a diverging 
length scale p"2HT8] . Such a “symmetric” scenario where 
£ diverges from both sides of the transition applies to the 
depinning transition m of an elastic manifold pushed 
through a disordered medium. Nevertheless, an alterna¬ 
tive scenario has been argued for in glassy systems with 
slowly-decaying interactions, predicting system-spanning 


avalanches (£ = oo) in the entire glass phase [EO]. Ap¬ 
plied to amorphous solids, this view suggests criticality 
for all stresses E < X max where plasticity occurs. This 
approach however lacks empirical support, and its conse¬ 
quences on failure near E max have not been investigated. 

In this Letter we show that as the stress is adiabat- 
ically increased in the solid phase, leading to a plas¬ 
tic strain e(X), the mean avalanche size ( S ) follows 
(S) ~ /{dY/de), where N is the system size and 

0 is an exponent that characterizes the stability of the 
structure m- This result confirms that avalanches are 
system-spanning (£ = oc) for all E < E max , and further 
implies an additional singularity as failure is approached, 
since dY/de 0 when E —)> X max . We suggest that 
data analysis used in the literature can mistakenly in¬ 
terpret this singularity as a diverging length scale. We 
also derive a scaling relation between 0 and exponents 
characterizing the statistics of avalanches. We test these 
predictions using elasto-plastic models [22[ |23], and show 
that they hold independently of the system preparation 
and of the presence of shear bands near E max , thus imply¬ 
ing that macroscopic flow localization and singularities in 
avalanche size are unrelated. 

Elastoplastic viewpoint: following [22H24] we model 
amorphous solids as consisting of N blocks, each char¬ 
acterized by a scalar local stress cq and a local failure 
threshold af 1 . The overall shear stress is E = Y li&i/N . 
Stability of i is achieved if || < erf 1 . Otherwise, the 
block is unstable: a plastic strain of magnitude Ae^ oc¬ 
curs on some time scale r c , leading to an overall incre¬ 
ment of plastic strain Ae = A ei/N. Such a plastic event 
also reduces stress locally by some amount A m /iA 
where /a is the elastic modulus, and affects stress in other 
locations via a long-range Eshelby field 5aj = G(rij)Aai 
[25] . which can in turn trigger new instabilities. For 
our numerics below, we choose the specific model de¬ 
scribed in [13 in two dimensions. Blocks then form a 
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bi-periodic square lattice, and the elastic propagator fol¬ 
lows approximatively Q{fij) oc cos(4</>)/r 2 where is the 
angle between the shear direction and r^. We choose 
<7 th = r c = i± = 1, and Acq = —cq + S, where 5 is a ran¬ 
dom number, uniformly distributed in [—0.1,0.1]. For 
these choices, E c ~ 0.53, and stability is easily expressed 
in terms of the variable aq = crj h — , and corresponds 

to Xi G [0, 2]. 

Such automaton models can be used to study the tran¬ 
sient regime toward failure. In what follows we use two 
quasi-static protocols. In the stress-control protocol, E 
is increased just sufficiently to trigger a single instabil¬ 
ity. E is fixed during the resulting avalanche, and is in¬ 
creased again only when this chain of events has stopped. 
The strain-control protocol is identical, except that E 
decreases during avalanches, proportionally to the plas¬ 
tic strain. Stress v.s. plastic strain curves for these two 
protocols are shown in Figjl] (from which the stress v.s. 
total strain 7 curves are easily deduced using the rela¬ 
tion Ay = Ae + AE/g). They essentially track each 
other macroscopically (although they differ microscopi¬ 
cally, see insets (c) and (d)) except when E reaches E max , 
if E max > E c . 

The transient qualitatively depends on the initial sta¬ 
bility of the system, characterized by the initial distribu¬ 
tion of local stability Po(x). If Pq{x) is narrow and de¬ 
pleted near x = 0 (corresponding to a very stable initial 
condition), transient shear bands occur; otherwise flow 
can remain homogeneous [22]. In Fig[l]we confirm these 
results using a broad and a narrow distribution Pq(x) 
(see S.I. for details). We further find that transient shear 
bands tend to occur if the stress-strain curve overshoots 
(although we did not investigate this correlation system¬ 
atically), as is sometimes reported [3], 27, 2S: and argued 
for in 0 [29 . In what follows we focus on avalanche-type 
response, for E below and approaching E max . 

Distribution of local distance to yield stress: mean-field 
models mug reveal that the distribution of local stabil¬ 
ity P(x) vanishes near x = 0 in a quasi-static shear at E c . 
In [21] some of us showed that stability indeed requires 
the presence of a pseudo-gap , i.e. P(x) ~ x e with 0 > 0, 
otherwise, any plastic event would eventually trigger an 
extensive rearrangement, and this argument also holds 
in the transient regime. 0 was measured in elasto-plastic 
models m and indirectly in MD simulations ED [32] 
both at E c and after a quench at E = 0, leading to con¬ 
sistent results. In Figj2]we extend these results to the 
transient regime. We find that 0 > 0 as predicted in [ 21 ]. 
However, the value of 6 turns out to be function of the rel¬ 
ative stress E/E max , while it converges to a well-defined 
value for large system size as shown in S.I. After some 
initial decay at very small E (not shown), the value of 0 
increases from 6 = 0.174 db 0.004 at E/E macc ~ 0.49 to 
the value 9 = 0.6 =b 0.004 at E = E max . This measure is 
consistent with the exponent obtained in the stationary 
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FIG. 1. Stress E v.s. plastic strain e curves for both 
strain (blue) and stress (red) controlled protocols, for (a) 
a broad initial distribution Pq(x) and (b) a narrow Po(x). 
In (a) E max = E c , whereas in (b) the stress overshoots and 
Emax > E c . Insets: spatial maps of plastic strain measured 
at different strain level. Highly (weakly) plastic regions are 
indicated in yellow (blue). Macroscopic shear localization oc¬ 
cur in (b) but not in (a). (c,d): Zooming in on the stress v.s. 
plastic strain curves, one observes microscopic differences be¬ 
tween the two protocols. 
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FIG. 2. (a) Distribution of local stability P{x) for E/E macc = 
0.49 to 0.99 in the stress-control protocol for N — 1024 2 in 
the case E max = E c . The dashed lines are direct fits of the 
form P(x) ~ x e , from which we extract 9. This quantity is 
shown in (b). The dashed line is the interpolation using a 
third order polynomial fit. 

The value of 9( E c ) was argued to control rheological 
properties in the flowing phase as E —»• E c from above|T3] 
and to imply system spanning avalanches for E < E maa , 
ED] . We now extend this latter argument to include the 
case where E —>> E maa , from below. 

Extreme value statistics implies that if P{x) ~ x e and 
the variables aq are independent, the least stable block 
must be at a distance x m i n ~ N-0+ 1 of an instability. 
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By definition, x min is the increment of stress that can be 
added before a new avalanche starts: the length of the 
vertical lines in Figjljd). Hence, following a finite stress 
increment of AH, a number M ~ AH/r m i n ~ AH7V0 + 1 
of avalanches are triggered. In elastoplastic models the 
avalanche size S is defined as the number of plastic 
events, which is approximately related to the total strain 
of the single avalanche Se , by S ~ NSe. Thus the total 
strain increase Ae must follow Ae = M(Se) = M(S)/N , 
where (S) is the mean avalanche size at stress H. We 
thus get: 


NAe _ AT®/( 1 + (? ) Ae v N e ^ 1+9 ") 
' ' “ M ~ AE > <9S/<9e 


( 1 ) 


where SH /de is the local slope of the stress-plastic strain 
curve, and the limit corresponds to AH —0. This cen¬ 
tral result indicates that (i) if H is increased in the solid 
phase, avalanches are system-spanning (£ = oo) even for 
H < H max , since their size is N dependent. Thus the 
system remains critical in the whole range H < H maa; 
as long as plastic flow occurs, i.e. 3H/3e < oo. (ii) 
Avalanches become larger as H ^ H max , as observed [ 8 ], 
since SH/de 0 at H max . 

Further scaling relations can be derived for the statisti¬ 
cal properties of transient avalanches for H < H max . We 
make the assumption that the distribution of avalanches 
P(S) is homogeneous, i.e. P(S) = S~ T f(S/ 5 C ), where 
the cut-off size scales as S c ~ L d f . Here df is the fractal 
dimension of avalanches, L is the linear system size, and 
N = L d , where d is the spatial dimension. From this 
distribution it is straightforward to compute the mean 
(S) ~ L d f( 2 ~ T \ Comparing with Eq.|l]), we get: 


d 6 

d~ f o TT 


( 2 ) 


A similar relation holds for stationary flow [13] , although 
in the transient regime exponents appear to depend con¬ 
tinuously on H. 

Finally, we introduce an exponent 7 defined as 
dS/de ~ (E max - E ) 7 for E close to S mas . Eq.Q then 
implies the scaling relation: 


(5) ~ (E max - E)- 7 a4t (3) 

These predictions are tested in Figj3] The inset of 
panel (a) shows that the mean avalanche size, as a func¬ 
tion of A = , grows with the system size even 

far from failure. The entire solid phase is critical, as ex¬ 
pected from Eq. 0 . Note that to test this equation, one 
must consider the fact that 0 = 0(H). In this figure we 
use for 0(H) the third order polynomial fit of Fig{2|b). 
Using these values for 0 a beautiful collapse is observed. 

The presence of system sized avalanches far from 
threshold has to be distinguished from the divergence 


observed close to the yield stress, (S) ~ (H max — H ) -7 at 
fixed N as implied by Eq.([3|. Figj3^a) is consistent with 
this relation and yields 7 « 1.1. According to its defini¬ 
tion, 7 can also be directly measured from the local slope 
of stress-strain curves, as is done in Figj3^b) where 7 ~ 1 
is found, consistent with Figj3ja). 7 = 1 means that the 
stress tends to H max exponentially fast. As shown in S.I, 
this appears to be valid also if the stress overshoots and 

^max ^ £ c . 

In Figj^c ,d), we measure df at A = (H max — 
S)/H max = 0.2, where 0 ~ 0.33, by collapsing the 
probability distribution of avalanche sizes, P(S) ~ 
S~ T f(S/S c ) with S c - L d f. We find d f « 0.77 and 
r « 1.35 in the stress-control case. Again, these val¬ 
ues perfectly agree with Eq. 0 . This result holds also 
for the strain-control protocol where we find scale-free 
avalanches with the same r and a similar fractal dimen¬ 
sion, df « 0.71. 



FIG. 3. (a) Collapse of the mean avalanche size as function of 
the proximity to failure A = E g ax ~ s , using the value of 0(H) 
extracted from Fig[2jb). The slope of the dotted line is —1.1. 
Here H max depends on the system size. The inset is the same 
quantity with no rescaling, (b) Local slope d'E/de vs H max — 
H, supporting 7 « 1, corresponding to dH/de ~ (H max — 
H) asymptotically. (c)&(d) Collapse of the distribution of 
avalanche size at a specific stress value corresponding to A = 
0.2 for the stress-control case(c) and strain-control case(d). 
We get df ~ 0.77(0.71) in the stress(strain)-control case, and 
r « 1.35 in both cases. All numerics are for d — 2 . 

Length scale: To further evidence the presence of a di¬ 
verging length scale throughout the solid phase, we study 
the strain map generated by a single avalanche, and con¬ 
sider the M = (S — 1)5/2 distances | R4 — Rj\ between 
all the blocks involved in the avalanche. We compute 
the distribution of these distances, and define P(R) as 
the average of these distributions among avalanches oc- 
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FIG. 4. (a) Distribution of avalanche extension P(R) for the stress-control protocol at A = 0.3. Collapse occur by rescaling 
distances with L, supporting that £ = L. Inset: direct measurement of £, defined as £ ~ (R). (b) C^R) at A = 0.3 for different 
Ae as indicated in legend, suggesting a length scale that depends on Ae. Inset: C 2 (R) for which no such dependence appears, 
(c) C 2 (R) computed for Ae = 10 -3 and varying A as indicated in legend, suggesting an increasing length scale as A —)> 0. 
Inset: C 2 (R) shows no such effect. 


curring at the same stress value in different samples (we 
choose to weight each avalanche by M in this average). 
We focus on avalanches occurring at a finite distance from 
failure, with A « 0.3. Assuming homogeneity we expect 
P(R) = (f)- I n FigJ^a), we confirm such a form, 

specifically: 

P(R) = (f) (4) 

We observe a similar scaling form in the strain-control 
simulation. These results confirm that £ ~ L, as further 
supported by the observation that (R) ~ L shown in 
inset. 

Our results are at odd with the conclusions of m, 
which report an increasing length scale in a stress-control 
simulation of granular media. We now suggest that their 
data may in fact be consistent with our views. In EH, a 
length scale is extracted by considering the fluctuations 
of the strain field obtained during some strain interval 
Ae, for different stress values E. This is a priori different 
from our analysis above which considers avalanches in¬ 
dividually. To clarify this point, we perform an analysis 
closer to theirs, where finite intervals of strain are con¬ 
sidered. We define a pair density function [33, C^R) as 
the probability that two local plastic events among the 
M « NAe ones in this interval are at a distance smaller 
than R. FigQb) shows C^i?) for A = 0.3 and varying 
Ae as indicated in legend. At first sight, one may think 
that a length scale can be extracted from ), but 
that this length is Ae-dependent. We find that this de¬ 
pendence, however, can be cured by removing the effect 
of the mean strain in our definition of C 2 (R ). We define 
0 2 {R) = C 2 {R) — C 2 {L)R 2 /L 2 , which is zero if the plastic 
events are homogeneous in space. As shown in the inset 
of FigQb), the characteristic length in C 2 (R ) does not 
depend on Ae. In FigQc) we show a similar analysis as 
the proximity to failure A is varied. From C 2 (R) it would 
appear that a length scale grows as A —0. However, as 


shown in inset this is an artefact of this analysis, as C 2 (R) 
shows a constant length scale of order L, consistent with 
our prediction £ = oc. Experimental measurements of 
the anisotropic part of the strain field also support that 
the correlation length is always large and weakly depends 
on m- Our views could be further tested by perform¬ 
ing a similar analysis in similar experimental izm and 
numerical EE! data. 

Conclusion: Ref. [20] argues that glassy systems whose 
elementary excitations display sufficiently long-range 
interactions (including electron glass, mean-field spin 
glasses or spheres at random close packing) must dis¬ 
play criticality for an entire range of fields or shear 
stress. This view has not yet been established experi¬ 
mentally. Our work supports that it holds in amorphous 
solids and granular materials, where it should be testable. 
Slowly-sheared granular material experiments have re¬ 
vealed avalanches with power-law statistics, but currently 
these studies have been limited to stationary flow [34] [35] 
(which, in addition to miss the transient behavior, may 
lead to additional complexity for granular materials due 
to the emergence of isostaticity [361 37i, a property how¬ 
ever absent in the transient [38]). Note that for stresses 
far from E max , a large system may be required to test 
our views, since one must have (S) 1 for Eq. 0 to 

hold. Our predictions may also apply in disordered crys¬ 
tals, where £ max however is not well-defined presumably 
due to work hardening [39| j40j . In [41] , the authors ob¬ 
serve numerically scale-free avalanches with (S) rsj N oa 
for a range of stresses. In our views that corresponds to 
0 « 0.67, a prediction that could be tested by measuring 
how the characteristic interval of stress with no plasticity 
vanishes with N. Finally, a central question for the future 
is what governs the value of the exponent #, which affects 
plasticity but also macroscopic rheological properties. 

Acknowledgments: It is a pleasure to thank E. 
DeGiuli, E. Lerner, B. Metzger, M. Muller, J. Weiss, 
L. Yan and S. Zapperi for discussions related to this 



























5 


work. MW acknowledges support from NSF CBET 
Grant 1236378 and MRSEC Program of the NSF DMR- 
0820341 for partial funding. 


[11 M. L. Manning, J. S. Langer, and J. M. Carlson, Physical 
review E 76 (2007). 

[2] B. Andreotti, Y. Forterre, and O. Pouliquen, Granular 
media: between fluid and solid (Cambridge University 
Press, 2013). 

[3] E. A. Jagla, Phys. Rev. E 76, 046119 (2007) 

[4] T. Divoux, M. A. Fardin, S. Manneville, and S. Lerouge, 
arXiv preprint arXiv: 1503.04130 (2015). 

[5] R. L. Moorcroft and S. M. Fielding, Physical review let¬ 
ters 110, 086001 (2013). 

[ 6 ] D. A. Lockner and N. M. Beeler, International Geo¬ 
physics 81, 505 (2002). 

[7] A. Amon, V. B. Nguyen, A. Bruand, J. Crassous, and 
E. Clement, Phys. Rev. Lett. 108, 135502 (2012) 

[ 8 ] A. Le Bouil, A. Amon, J.-C. Sangleboeuf, H. Orain, 
P. Besuelle, G. Viggiani, P. Chasle, and J. Crassous, 
Granular Matter 16, 1 (2014). 

[9] A. Le Bouil, A. Amon, S. McNamara, and J. Crassous, 
Physical review letters 112 , 246001 (2014). 

[10] F. Gimbert, D. Amitrano, and J. Weiss, EPL (Euro¬ 
physics Letters) 104, 46001 (2013). 

[ 11 ] A. Argon, Acta Metallurgica 27, 47 (1979) 

[12] O. Pouliquen, Physical review letters 93, 248001 (2004). 

[13] J. Lin, E. Lerner, A. Rosso, and M. Wyart, Proceedings 
of the National Academy of Sciences 111, 14382 (2014). 

[14] P. Olsson and S. Teitel, Phys. Rev. Lett. 99, 178001 
(2007). 

[15] K. Martens, L. Bocquet, and J.-L. Barrat, Phys. Rev. 
Lett. 106, 156001 (2011) 

[16] G. During, E. Lerner, and M. Wyart, Physical Review 
E 89, 022305 (2014). 

[17] K. M. Salerno, C. E. Maloney, and M. O. Robbins, Phys. 
Rev. Lett. 109, 105703 (2012) 

[18] A. Lemaitre and C. Caroli, Phys. Rev. Lett. 103, 065501 
(2009) 

[19] D. S. Fisher, Physics Reports 301, 113 (1998) 

[20] M. Muller and M. Wyart, Annual Review of Condensed 
Matter Physics 6 (2015). 

[21] J. Lin, A. Saade, E. Lerner, A. Rosso, and M. Wyart, 
EPL (Europhysics Letters) 105, 26003 (2014). 

[22] G. Picard, A. Ajdari, F. Lequeux, and L. Bocquet, Phys¬ 
ical Review E 71, 010501 (2005). 

[23] J.-C. Baret, D. Vandembroucq, and S. Roux, Phys. Rev. 
Lett. 89, 195506 (2002) 

[24] P. Hebraud and F. Lequeux, Phys. Rev. Lett. 81, 2934 
(1998), 

[25] G. Picard, A. Ajdari, F. Lequeux, and L. Bocquet, The 
European Physical Journal E 15, 371 (2004) 

[26] D. Vandembroucq and S. Roux, Physical Review B 84, 
134210 (2011). 

[27] Y. Shi, M. B. Katz, H. Li, and M. L. Falk, Physical 
review letters 98, 185505 (2007). 

[28] R. L. Moorcroft, M. E. Cates, and S. M. Fielding, Phys¬ 
ical review letters 106, 055502 (2011). 

[29] S. M. Fielding, Reports on Progress in Physics 77, 102601 
(2014). 


[30] A. Lemaitre and C. Caroli, arXiv preprint 
arXiv:0705.3122 (2007). 

[31] S. Karmakar, E. Lerner, and I. Procaccia, Phys. Rev. E 
82,055103 (2010) 

[32] K. M. Salerno and M. O. Robbins, Physical Review E 
88 , 062206 (2013). 

[33] L. Girard, D. Amitrano, and J. Weiss, Journal of Statis¬ 
tical Mechanics: Theory and Experiment 2010, P01013 
( 2010 ). 

[34] N. W. Hayman, L. Ducloue, K. L. Foco, and K. E. 
Daniels, Pure and applied geophysics 168, 2239 (2011). 

[35] M. Bretz, R. Zaretzki, S. B. Field, N. Mitarai, and 
F. Nori, EPL (Europhysics Letters) 74, 1116 (2006). 

[36] E. DeGiuli, G. During, E. Lerner, and M. Wyart, Phys. 
Rev. E 91, 062206 (2015) 

[37] E. Lerner, G. During, and M. Wyart, Proceedings of the 
National Academy of Sciences 109, 4798 (2012). 

[38] N. P. Kruyt, Comptes Rendus Mecanique 338, 596 

( 2010 ). 

[39] M.-C. Miguel, A. Vespignani, S. Zapperi, J. Weiss, and 
J.-R. Grasso, Nature 410, 667 (2001) 

[40] F. F. Csikor, C. Motz, D. Weygand, M. Zaiser, and 
S. Zapperi, Science 318, 251 (2007) 

[41] P. D. Ispanovity, L. Laurson, M. Zaiser, I. Groma, S. Zap¬ 
peri, and M. J. Alava, Phys. Rev. Lett. 112 , 235501 
(2014). 



6 


APPENDIX 


(a) Dependence on initial conditions 


Different transient behaviors can be generated by con¬ 
sidering different initial conditions for the distribution 
of local stability Po(x), which we build as follows. We 
set initially Xi = 1 — JT G(rij)Axj, where A Xj is a ran¬ 
dom variable uniformly distributed in [0, Ax]. Such ini¬ 
tial conditions ensure that each line and column present 
the same mean x, which reduces finite size effects m- 
The dynamics is then run at fixed stress, until all sites 
are stabilized. We choose Ax = 3.0 (broad), and Ax = 
0.7(narrow), as for Fig.l in the main text. The two cor¬ 
responding distributions are shown in Fig |A.l[ a). 



FIG. A.l. (a) The initial distributions Po(x) before shear 
starts. Blue solid line corresponds to Ax = 3.0 (broad, no 
overshoot) whereas green dashed line corresponds to Ax = 
0.7 (narrow, overshoot), (b) Rescaled mean avalanche size 
(. S)/N 0-2 for the overshoot case(Ax = 0.9), where (S) is again 
observed to be size -dependent, and 7 ~ 1. The dased line 
has a slope — 1 . 


Our predictions, including the result of Eq.(l), appear 
to hold true independently of the system preparation, as 
shown in Fig |A.l[ b) testing Eq.(l) in the overshoot case. 
Although we did not perform a detailed measurement of 
the exponent 0(E) in that case, we find that (S) is indeed 
system size dependent, and again find 7 ~ 1. 


(b) Finite size effects on P(x) and 0(E) 

Fig |A.2| shows P(x) for different system sizes and 
four different stress values, supporting that the value 
we report for 0(E) in the main text, Fig.2(b), indeed 
approximates the true value in the thermodynamic limit. 



FIG. A.2. P(x) as a function of E and system size in the 
stress-controlled protocol. 













